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Abstract 



Measurements of jet production rates in association with W and Z bosons for jet trans- 
verse momenta above 30 GeV are reported, using a sample of proton-proton collision 
events recorded by CMS at ^/s = 7 TeV, corresponding to an integrated luminosity 
of 36pb^^. The study includes the measurement of the normalized inclusive rates of 
jets cr(y+ > n jets)/t7'(V), where V represents either a W or a Z. In addition, the 
ratio of W to Z cross sections and the W charge asymmetry as a function of the num- 
ber of associated jets are measured. A test of Berends-Giele scaling at v^s = 7 TeV 
is also presented. The measurements provide a stringent test of perturbative-QCD 
calculations and are sensitive to the possible presence of new physics. The results are 
in agreement with the predictions of a simulation that uses explicit matrix element 
calculations for final states with jets. 
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1 Introduction 

The study of jet production in association with a W or Z vector boson (denoted by 'V in this 
paper) provides a stringent test of perturbative QCD calculations. The presence of a vector 
boson provides a clear signature of the process and allows comparison of different scattering 
amplitudes with respect to inclusive multijet production, which is dominated by gluon scatter- 
ing. In addition, a precise measurement of the W (Z) + n jets cross section is essential since the 
production of vector bosons with jets constitutes a background in searches for new physics and 
for studies of the top quark. At present, next-to-leading-order (NLO) predictions are available 
for V + n jets, with n up to four [1-5J. 

The W+ to cross-section ratio for an associated jet multiplicity is also predicted by the 
theory [6J and larger than unity. Deviations of this ratio from the expected value may point to 
new processes that produce W bosons in the final state. 

The CDF, DO and ATLAS Collaborations have reported measurements of W (Z) + n jets pro- 
duction in llTHTTll. 

This paper presents results obtained with the 2010 data sample of the Compact Muon Solenoid 
(CMS) experiment at the Large Hadron Collider (LHC), based on an integrated luminosity of 
35.9 ± 1.4 pb^^ [12J, collected in proton-proton collisions at 7 TeV. To reduce systematic uncer- 
tainties associated with the integrated luminosity measurement, the jet energy scale (JES), the 
lepton reconstruction, and the trigger efficiencies, we measure the V + n jets cross sections rel- 
ative to the inclusive W and Z cross sections, cr(V+ > n jets) / c(V). A clear advantage of using 
ratios of cross sections is also that theoretical uncertainties tend to cancel, improving the ro- 
bustness of the results. We also measure the ratio of the W to Z cross sections and the W charge 
asymmetry as a function of the associated jet multiplicity. Finally, we test the Berends-Giele 
scaling hypothesis |[13| . To compare the measurements with the predictions of the theoretical 
calculations, the results are presented at "particle level", unfolding the detector efficiency and 
resolution. 

2 The CMS Experiment 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal di- 
ameter, providing an axial magnetic field of 3.8 T. Within the field volume are a silicon pixel 
and strip tracker, an electromagnetic calorimeter (ECAL), and a brass/ scintillator hadron cal- 
orimeter (HCAL). Muons are detected in gas-ionization detectors embedded in the steel return 
yoke. 

A right-handed coordinate system is used in CMS, with the origin at the nominal interaction 
point, the x-axis pointing to the centre of the LHC ring, the y-axis pointing up (perpendicular 
to the LHC plane), and the z-axis along the anticlockwise-beam direction. The polar angle 9 
is measured from the positive z-axis and the azimuthal angle (p is measured from the positive 
X-axis in the x-y plane. The pseudorapidity is defined asrj = — ln[tan(0/2)]. 

The inner tracker contains 1440 silicon pixel and 15148 silicon strip detector modules. It mea- 
sures charged particle trajectories in the pseudorapidity range \rj\ < 2.5. The ECAL consists 
of nearly 76000 lead tungstate crystals that provide coverage for \rj\ < 1.479 in a cylindrical 
barrel region (EB) and 1.479 < \rj\ < 3.0 in the two endcap regions (EE). Preshower detectors, 
each consisting of two planes of silicon sensors interleaved with a total of 3 radiation length 
of lead, are located in front of both EEs. The HCAL is a sampling device with brass as the 
passive material and scintillator as the active material. The combined calorimeter cells are 
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3 Data and Simulation Samples 



grouped in projective towers of granularity A/y x = 0.087 x 0.087 at central rapidities and 
/^rj X A^ « 0.17 X 0.17 at forward rapidities. Muons are detected in the pseudorapidity range 
\rj\ < 2.4, with detection planes based on three technologies: drift tubes, cathode strip cham- 
bers, and resistive plate chambers. In addition to the barrel and endcap detectors, CMS has 
extensive forward calorimetry. The first level (LI) of the trigger system, composed of custom 
hardware processors, is designed to select event candidates in less than 3.2 }is using informa- 
tion from the calorimeters and muon detectors. The high level trigger (HLT) processor farm 
further reduces the event rate to a few hundred hertz, before data storage. A more detailed 
description of CMS can be found elsewhere |(T4). 

3 Data and Simulation Samples 

The LI trigger system selected electrons with an energy deposit in the ECAL of at least 5 GeV 
or 8 GeV, depending on the luminosity conditions, and muons with a transverse momentum 
exceeding 7 GeV. The events were then filtered by the HLT with algorithms that evolved in 
response to the rapid rise of the LHC luminosity during 2010. The px threshold of the electrons 
and muons was adjusted periodically, to cope with the increasing instantaneous luminosity. 
The largest sample of electrons was collected with an online requirement > \7 GeV. For 
muons, most data were collected with a threshold of px ^ 15 GeV. 

Simulated Monte Carlo (MC) samples are used for data /simulation comparison and to unfold 
the jet multiplicity distributions. Backgrounds are estimated from data, as explained below. 

Samples of simulated events with a W or a Z boson are generated with the MadGraph 4.4.13 
IflSJ event generator, interfaced to the PYTHIA 6.422 [16| program for parton shower simula- 
tion. The set of parton distribution functions used is CTEQ6L1 |[T7|. The MadGraph gen- 
erator produces parton-level events with a vector boson and up to four partons on the basis 
of a matrix-element calculation. This sample serves as the baseline for comparisons with data. 
Additional samples of W and Z events are generated with pythi A. For jet multiplicities greater 
than one, MadGraph is expected to be more accurate since it uses the exact matrix-element 
calculation, while in PYTHIA only the hardest emission reproduces the exact matrix-element 
calculation. Top-pair (tf) and single-top production processes are generated with MadGraph. 
Processes of multijet, 7+jets, b-hadron and c-hadron decays to final states with electrons and 
muons are generated with PYTHIA. 

The full list of simulated samples is given in Table [ij In order to compare with the data dis- 
tributions, the simulation samples are normalized to the cross sections times the integrated 
luminosity, using or NLO or next-to-next-to-leading-order (NNLO) cross sections IIT8] - t20l , or 
the leading-order (LO) cross sections from the MC generator, as reported in the table. 

Generated events are processed through a full detector simulation program based on GEANT4 
pn [22| , followed by a detailed emulation of the trigger and event reconstruction. Minimum- 
bias events are superimposed on the generated events to reproduce the distribution of multiple 
proton-proton collisions in the same bunch crossing (pileup) observed in 2010. A signal sample 
without pileup is used for purpose of comparison. The PYTHIA parameters for the underlying 
event are set to the Z2 tune [23], which is a modification of the Zl tune described in Ref. [|24|. 
Comparisons are also made to the D6T tune p5|. 
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Table 1: Summary of simulated datasets for the various signal and background processes used 
in this analysis. The requirements applied to leptons pj and t], dilepton invariant mass {Mu), 
and transverse momentum of the hard interaction (pV), are shown, and the corresponding cross 
section is given. 



Process 


Generator 


Kinematic selection 


a(ph) 




MadGraph 


no selection 


3.1x10" (NNLO) 


z ^ i+i- 


MadGraph 


Mu > 50 GeV 


3.0 xlO^ (NNLO) 


tt 


MadGraph 


no selection 


1.6x10^ (NLO) 


Single-top tW channel 


MadGraph 


no selection 


1.1 xlO^ (LO) 


Single-top s and t channels 


MadGraph 


no selection 


3.5 (NLO) 


W ^> ei/ 


PYTHIA 


\%'\ <2.7 


8.2 xlO^ (NNLO) 


W ^> /iy 


PYTHIA 


h'\ <2.5 


7.7 xlO^ (NNLO) 


W TV 


PYTHIA 


no selection 


1.0 X 10* (NNLO) 


Z ^ 


PYTHIA 


Mu > 20 GeV 


5.0 xlO^ (NNLO) 


Inclusive }i QCD multijet 


PYTHIA 


Pt > 20 GeV, pl^ > 10 GeV, \t],,\ < 2.5 


3.4 xlO^ (LO) 


EM-enriched QCD multijet 


PYTHIA 


20 GeV < pr < 170 GeV 


5.4x10^ (LO) 


b/c e 


PYTHIA 


20 GeV < pr < 170 GeV 


2.6 xlO^ (LO) 


7+jet 


PYTHIA 


no selection 


8.5x10^ (LO) 



4 Signal Selection 

Signal selection begins with the identification of a charged lepton, either an electron or a muon, 
with Pj > 20 GeV. This lepton, which will be called the "leading lepton", must geometrically 
match the object that triggered the event readout. 

For electron candidates, we require that the ECAL cluster lies in the fiducial region < 2.5, 
with the exclusion of the region 1.4442 < \f]\ < 1.566. This exclusion allows us to reject elec- 
trons close to the barrel/ endcap transition and electrons in the first endcap trigger tower, which 
lies in the shadow of cables and services. A series of quality requirements, following the stan- 
dard established by the measurement of the inclusive W and Z cross sections 1 26 1 , are then 
applied to the electron, as briefly described below. 

For each electron candidate, a supercluster in ECAL is defined in order to correct for the po- 
tential underestimation of the energy due to bremsstrahlung. Thus the electron cluster is com- 
bined with the group of single clusters attributed to bremsstrahlung photons generated in the 
material of the tracker. Electrons are first selected based on the spatial matching between the 
ECAL supercluster and the silicon detector track in the ?/ and (p coordinates, on the supercluster 
energy distribution in the rj direction, and on the energy leakage into the HCAL detector. To 
reduce the contamination from converted photons, a minimal track transverse impact param- 
eter significance is required. Electrons are rejected if no associated track hits are found in the 
first tracker layers or if a conversion partner candidate is found. 

To reduce further the contamination from misidentified electrons and hadronic decays, we se- 
lect electrons isolated from hadronic activity. This selection is based on maximum values al- 
lowed for three isolation variables. The variables are computed relative to the electron Ej and 
consist of the sums of track pj in the tracker, energy deposits in ECAL, and energy deposits 
in HCAL, respectively. The sums are computed inside the cone AR = + (A(/>)2 < 0.3, 

where A// and Acp are the differences in the pseudorapidity and azimuthal angle between the 
electron and the track or energy deposit, with an inner exclusion region that removes the elec- 
tron contribution. 

For the leading electron, the values of the different quality requirements are chosen such that 
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5 Jet Rates 



they correspond to a lepton efficiency of about 80%, as evaluated with the MadGraph +PYTHIA 
simulated sample described in Section |6] 

After identifying the leading electron, we search for a second electron candidate called the 
"second leading electron", within the ECAL fiducial volume and with pj > 10 GeV. The quality 
requirements for the second leading electron are tuned to provide an electron efficiency of 
about 95%, as evaluated with the MadGraph +PYTHIA simulated sample. If such a second 
leading electron is found, and its invariant mass with the first leading electron Mfj lies between 
60 GeV and 120 GeV, the event is assigned to the Z + jets sample. If such a second leading 
electron is not found, the event is assigned to the W + jets sample, thereby ensuring that there 
is no overlap between the two samples. Events including a muon with pj > 15 GeV and 
|?7 1 < 2.4 are rejected from the W + jets electron sample to reduce tt contamination. 

The muon reconstruction and identification are identical to that used for the measurement of 
the W and Z cross sections [26 J. A relative isolation variable, I = J^ipj^"^ + E^^^^ + Ef^^^)/ p^ 
is defined, which includes the pj for tracks, ECAL, and HCAL towers in a cone AK < 0.3 
around the muon direction. The muon and its energy deposits are excluded from the sum. 
A muon is considered to be isolated if / < 0.15. The V + jets muon event selection starts by 
requiring the presence of an isolated muon in the region \rj\ < 2.1 with pi > 20 GeV. It must 
be a high-quality muon as described in [26| with an impact parameter in the transverse plane 
\dxy\ < 2 mm, to suppress cosmic-ray background. As for the electron channel, we then search 
for a second leading muon with pj > 10 GeV and \ri\ < 2.4, such that the dimuon invariant 
mass lies between 60 GeV to 120 GeV. If such a second leading muon is (is not) found, the event 
is assigned to the Z + jets (W + jets) sample. An electron veto is not applied in the selection of 
W decays into muons as it would significantly lower the efficiency and increase the systematic 
uncertainty, because of fake electrons. 

For the W + jets samples in both decay modes, we compute the missing transverse energy £t 
using a particle-flow (PF) algorithm pTil that reconstructs individually each particle in the event 
based on information from all relevant subdetectors. We then use it to calculate the transverse 
mass, Mj = -^2^x^1(1 ~ cos A^), where Atp is the angle in the x-y plane between the directions 
of the lepton pj and the ^j, and select events with Mj > 20 GeV. 

For the charge asymmetry measurement, the charge assignment of the lepton is used to select 
W+ + jets and + jets candidates. For electrons, three algorithms are used to determine 
their charge 128] : the curvature of the electron track reconstructed by a Gaussian-sum-filter 
algorithm II29II , the curvature of the silicon detector track associated with the electron, and 
the difference in (p between the track direction at the interaction vertex and the supercluster 
measured by the ECAL. In order to reduce uncertainties related to charge misassignment, we 
reject the event if the charge value differs amongst the three methods. 

5 Jet Rates 

Jets are reconstructed from the particle collection created with the particle-flow algorithm and 
are formed with the anti-A;^ clustering algorithm [30] with a size parameter of R = 0.5. Jet 
energy corrections (JEC) are applied to account for the jet energy response as a function of rj 
and Ex USB. 

We require the jet to have \}]\ < 2.4 to be in the tracker acceptance, and Ej > 30 GeV. Jets are 
required to satisfy identification criteria that eliminate jets originating from or seeded by noisy 
channels in the hadron calorimeter 1321 . 
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The pileup and the underlying event affect the jet counting by contributing additional energy to 
the measured jet energy and therefore "promoting" jets above the Ej threshold for jet counting. 
To minimize the uncertainty due to pileup and the underlying event, the Ej threshold is set at 
30 GeV. The average number of pileup events in the data sample is 2.7. Their effect is taken into 
account by evaluating event-by-event the energy not related to the hard-interaction activity 
||33l|34l. This amount is subtracted from each jet ||3T1| . Removing the average energy due to 
pileup does not remove the jets from pileup interactions. In a simulated sample, this additional 
contribution to the jet count, however, is found to be negligible for the chosen jet Ej threshold. 

Electrons can be reconstructed as jets or can overlap with a jet. Therefore, jets that fall within 
AR < 0.3 of an electron from W or Z decay are not included in the jet count. Muons can also 
overlap with a jet, thus, each selected muon is matched to one of the particles reconstructed by 
the PF algorithm and excluded from jet clustering. 

One of the most important backgrounds in the W sample at high jet multiplicity comes from 
tt events. These contain two b-quark jets. We count the number of b-tagged jets, n^J^^^^'^, 
with a tagging algorithm that requires at least two tracks in the jet with a significance on the 
transverse impact parameter greater than 3.3. The algorithm parameters chosen correspond to 
a working point with an efficiency of about 62% and a mistag rate of about 2.9%. The efficiency 
is measured from data using a sample of tt events with fully leptonic final states. The mistag 
rate is averaged over light jets in simulated W and top events, and corrected for measured 
differences between data and simulation Il35l . The nj^j^^^^*^ value is then used in the fitting 
method to extract the number of signal events, as described in Section [/[ 

The observed transverse energy distributions of the leading jet in the W+ > 1 jet and Z+ > 
1 jet samples are shown in Figs. [T] and |2| respectively. Figures |3] and |4] show the same distribu- 
tions for the second-leading jet in the W+ > 2 jets and Z+ > 2 jets samples. In addition to the 
selection described in Section|4| for the transverse energy distributions, W boson candidate are 
required to have Mj > 50 GeV, thus reducing significantly the QCD multijet contamination. 
The tt background is larger in the muon than in the electron charmel because an electron veto 
is not applied, as explained in Section|4j The data are in good agreement with the MadGraph 
+ PYTHIA parton shower predictions normalized to the NNLO cross sections. 

The selected events are assigned to bins of jet multiplicity by counting the number of jets with 
£t > 30 GeV. The number of events with a W (Z) candidate and at least one jet are 43561 (1648) 
in the electron channel, and 32496 (2339) in the muon channel. The observed distributions 
of the exclusive numbers of reconstructed jets in the W and Z samples are shown in Figs. |5] 
and |6j respectively. The distributions from simulation are also shown. They are normalized to 
the integrated luminosity using the cross-section values in Table [ij Comparisons between the 
data and simulation distributions show that there is a very good agreement, even prior to any 
attempt to fit the signal and background content of the data. For the W sample, QCD multijet 
processes constitute the most significant source of background in the lower jet-multiplicity bins, 
while top quark production dominate in the higher jet-multiplicity bins. For Z events with one 
or more jets, the top quark production dominates the background, which overall is very small. 



6 Acceptance and Efficiency 

In order to provide model-independent results, we do not correct for the detector acceptance, 
but rather quote the results within the acceptance, as defined by the lepton and jet fiducial and 
kinematic selections given above. 
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6 Acceptance and Efficiency 



CMS 



CMS 



36 pb ' at vs = 7 TeV" 
Mr > 50 GeV 

• data 
□ W^ev 

n top 

^1 other backgrounds 




36 pb^ at \fs = 7 TeV" 
Mr > 50 GeV 

• data 
I I W^|iv 
^1 top 

^1 other backgrounds 




100 150 200 250 
leading jet Ej [GeV] 



o 

■D 



2 ■ 

1 



.. •■•'■•■.■♦♦♦.♦♦■♦♦vtf^f^-^-f^i|^-f 



50 100 150 200 250 
leading jet Et [GeV] 



Figure 1: Distributions of the Ej for the leading jet in the W+ > 1 jet sample for the electron 
channel (left) and for the muon channel (right), before the requirement of £t > 30 GeV (shown 
by the vertical dotted line) is imposed for counting jets. Points with error bars are data, his- 
tograms represent simulation of the signal (yellow), tf and single top backgrounds (orange) and 
other backgrounds (purple). The last includes multijet events, Z events, and events in which 
the W decays to a channel other than the signal channel. Error bars represent the statistical 
uncertainty on the data. The ratio between the data and the simulation is shown in the lower 
plots. 
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Figure 2: Distributions of the Ej for the leading jet in the Z+ > 1 jet sample for the electron 
channel (left) and for the muon channel (right), before the requirement of Ej > 30 GeV (shown 
by the vertical dotted line) is imposed for counting jets. Points with error bars are data, his- 
tograms represent simulation of the signal (yellow), tt and single top backgrounds (orange) and 
other backgrounds (purple). The last includes multijet events, W events, and events in which 
the Z decays to a channel other than the signal channel. Error bars represent the statistical 
uncertainty on the data. The ratio between the data and the simulation is shown in the lower 
plots. 
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Figure 3: Distributions of the Ej for the second leading jet in the W+ > 2 jets sample for the 
electron channel (left) and for the muon channel (right), before the requirement of Ej > 30 GeV 
(shown by the vertical dotted line) is imposed for counting jets. Points with error bars are data, 
histograms represent simulation of the signal (yellow), tf and single top backgrounds (orange) 
and other backgrounds (purple). The last includes multijet events, Z events, and events in 
which the W decays to a channel other than the signal channel. Error bars represent the sta- 
tistical uncertainty on the data. The ratio between the data and the simulation is shown in the 
lower plots. 
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Figure 4: Distributions of the Ej for the second leading jet in the Z+ > 2 jets sample for the 
electron channel (left) and for the muon channel (right), before the requirement of Ej > 30 GeV 
(shown by the vertical dotted line) is imposed for counting jets. Points with error bars are 
data, histograms represent simulation of the signal (yellow), tt and single top backgrounds 
(orange) and other backgrounds (purple). The last includes multijet events, W events, and 
events in which the Z decays to a channel other than the signal channel. Error bars represent 
the statistical uncertainty on the data. The ratio between the data and the simulation is shown 
in the lower plots. 
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6 Acceptance and Efficiency 
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Figure 5: Exclusive number of reconstructed jets in events with W — > ei/ (left) and W — > f/v 
(right). Points with error bars are data, histograms represent simulation of the signal (yellow), 
tt and single top backgrounds (orange) and other backgrounds (purple). Error bars represent 
the statistical uncertainty on the data. These distributions have not been corrected for detector 
effects or selection efficiency. The ratio between the data and the simulation is shown in the 
lower plots. 
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Figure 6: Exclusive number of reconstructed jets in events with Z — )■ e+e^ (left) and Z 
(right). Points with error bars are data, histograms represent simulation of the signal (yellow), 
tt and single top backgrounds (orange) and other backgrounds (purple). Error bars represent 
the statistical uncertainty on the data. These distributions have not been corrected for detector 
effects or selection efficiency. The ratio between the data and the simulation is shown in the 
lower plots. 
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The efficiencies for lepton reconstruction, identification, isolation, and trigger are obtained 
from Z + jets data samples. The sample for the measurement of a given efficiency contains 
events selected with two lepton candidates. One lepton candidate, called the "tag", satisfies all 
selection requirements, including the matching to a trigger object. The other lepton candidate, 
called the "probe", is selected with criteria that depend on the efficiency being measured. The 
invariant mass of the tag and probe lepton candidates must fall in the range [60-120] GeV. The 
signal yields are obtained for two exclusive subsamples of events in which the probe lepton 
passes or fails the selection criteria considered. Fits are performed to the invariant-mass distri- 
butions of the "pass" and "fail" subsamples, including a term that accounts for the background. 
The measured efficiency is deduced from the relative levels of signal in the "pass" and "fail" 
subsamples. An estimate of the systematic uncertainty on the efficiency is obtained by varying 
the shape used to fit the Z signal among several functional forms tested on the simulation. The 
lepton selection efficiency is the product of three components: the reconstruction efficiency, the 
identification and isolation efficiency, and the trigger efficiency. Each of these efficiencies is cal- 
culated as a function of the jet multiplicity in the event. The efficiency of additional selection 
requirements applied to W events is computed from the simulation. 

For electrons we find that the efficiency is about 70% (60%) for the W + jets (Z + jets) sig- 
nal events with variations of a few percent across different jet multiplicity bins, which are 
accounted for in the measurement. The uncertainty ranges from about 1% up to about 5%, 
increasing with the number of jets in the event. The additional requirement that the three 
methods to measure the electron charge agree, which is applied for the charge asymmetry 
measurement, has an efficiency of 97%. 

For muons, the efficiencies are measured as a function of px, rj and the jet multiplicity. The 
efficiency dependence on pj and // is measured in events with zero and one jet. For events 
with more than one jet, where the sample size is insufficient for a (px/ v) dependent efficiency 
to be measured, an average efficiency is computed and the {pj, rj) dependence measured in 
the n = 1 bin is assumed. The isolation efficiency is found to have the largest dependence 
on multiplicity, while the trigger, reconstruction, and identification efficiencies are consistent 
with being constant. The average efficiency is 85% and 86% for W and Z events, respectively. 
They are similar because the efficiency for selecting the second muon from Z candidates and 
for passing the Mj criteria with W candidates are both near unity. While the second muon 
efficiency in Z events is fairly flat across different jet multiplicity bins, the efficiency of the Mj 
requirement for W candidates decreases to 90% for high jet multiplicity. The uncertainty on the 
muon efficiency ranges from about 3% up to about 15%, increasing with the number of jets in 
the event. 

Electron charge misidentification is estimated with simulated samples and data. It depends on 
the lepton rj, with the largest misidentification in the endcap regions. The probability ranges 
from 4 X 10^^ up to 3 x lO^'^ in the simulation, and from 1 x lO^'^ up to 4 x 10^^ in the data. 
The difference between simulation and data is taken as uncertainty. The muon charge misiden- 
tification has been found to be smaller than 10^^ in cosmic-ray data [36 1, therefore yielding a 
negligible systematic uncertainty. Using the tag-and-probe method, the ratio of the selection 
efficiency for electrons to positrons is measured to equal unity within the statistical precision 
of 1.4% [28J. The efficiency ratio between f/+ and is found equal to one within 2% [28J. 
Effects related to charge measurements are propagated as systematic uncertainties to the final 
W charge asymmetry results. 
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7 Signal Extraction 
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Figure 7: Dilepton mass for the Z + 1 jet samples, in the electron channel (left) and the muon 
channel (right). Points with error bars are data. The fit result for the signal is shown by the 
yellow-filled area. The background is too small to be visible in the figure. 

7 Signal Extraction 

The signal yield is estimated using an extended maximum-likelihood fit to the dilepton invari- 
ant mass (Mie) distribution for the Z + jets sample, and to the transverse mass (Mj) distribu- 
tion for the W + jets sample; the number of observed events is included in the likelihood fit 
as a constraint on the normalization. The probability distribution functions are asymmetric 
Gaussians with tails. Their parameters are derived from the simulation or, for the background, 
from control data samples with inverted identification (isolation) requirement on the electrons 
(muons). 

For the Z event sample, the contamination from the background processes, dominated by tt and 
W + jets, is small and does not produce a peak in the Mfj distributions, so the distribution 
is taken to be the sum of two components, one for the signal and one that accounts for all 
background processes. 

For the W sample, the background contributions can be divided into two components, one 
which exhibits a peaking structure in Mj, dominated by tt, and another which does not, dom- 
inated by QCD multijet events. We perform a two-dimensional fit to the Mj distribution and 
the number of b-tagged jets, n^J^^^'^'^. The Mj distribution allows the statistical separation of 

the signal from the non-peaking backgrounds, while nj^^'^^^'^'^ distinguishes the signal and the 
background from tt. The likelihood function is built under the assumption that there are no 
b jets in the signal events. This implies that a fraction of W events produced in association 
with heavy-flavour jets, i.e. the fraction with at least one heavy-flavour jet in the acceptance, is 
counted as background. Considering the statistical precision of the measurement, this assump- 
tion has a negligible effect on the W + jets cross-section result. 

The fits are done using the jet multiplicity bins for n < 3; in contrast, the jet counting is done 
inclusively for the last bin of jet multiplicity, i.e. n > 4. Examples of fits for Z + 1 jet are shown 
in Fig.jzj Figuresjs and|9jshow fits in Mj and nj^^'^^^^*^ projections for the W + n jets (n=l and 
n=3) channel. The presence of the top quark background is evident from comparing the n = 1 
and n = 3 multiplicity bins. The fit is repeated separately on the W+ + jets and + jets 
samples for the charge-asymmetry measurement. 

In the electron channel, the observed exclusive V + jets yields determined from the fit are cor- 
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Figure 8: Fit results for the W(ei/) + n jets sample with n = 1 (upper row) and n = 3 (lower 
row). For each row, the Mt projection is shown on the left, while the nj^j'^^^'^'^ projection is 
shown on the right. Points with error bars are data. Fit results are shown by the colour-filled 
areas for the signal process (yellow), non-top backgrounds (purple), and top backgrounds (or- 
ange). 
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Figure 9: Fit results for the W(/^i/) + n jets sample with n = 1 (upper row) and n = 3 (lower 
row). For each row, the Mt projection is shown on the left, while the nj^j'^^^'^'^ projection is 
shown on the right. Points with error bars are data. Fit results are shown by the colour-filled 
areas for the signal process (yellow), non-top backgrounds (purple), and top backgrounds (or- 
ange). 
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rected a posteriori for electron efficiencies, which are discussed in Section |6] In the muon chan- 
nel, efficiencies are available in bins of the lepton pj, tj and the jet multiplicity. This allows an 
efficiency-weighted fit to be performed, which returns efficiency-corrected yields. 

8 Unfolding 

In order to estimate the scaling behaviour of the jets at the particle-jet level, we apply an unfold- 
ing procedure that removes the effects of jet energy resolution and reconstruction efficiency. 
The migration matrix, which relates a number n' of produced jets at the particle level to an 
observed number n of reconstructed jets, is derived from simulated samples of Z + jets and 
W + jets. The unfolding procedure takes into account the statistical and systematic uncertain- 
ties in the signal yields. 

The migration matrices are derived from simulated events with leptons and jets within the px 
and Tj acceptance of the analysis. The simulation is performed using the MadGraph generator 
with the Z2 tune and incorporates event pileup. 

Two unfolding methods are employed. The baseline method is the "singular value decompo- 
sition" (SVD) algorithm [37]. As a crosscheck, the iterative method [38J is also applied. Both 
algorithms require a regularization parameter to prevent the statistical fluctuations in the data 
from appearing as structure in the unfolded distribution. For the SVD method, the regulariza- 
tion parameter is chosen to be fcsvD = 5, corresponding to the number of bins. This algorithm 
corresponds to an inversion of the migration matrix, and gives the most accurate uncertainty 
estimate. For the iterative algorithm, the regularization parameter /cBayes = 4 is used, as sug- 
gested in Ref. |i39j|. 

9 Systematic Uncertainties 

One of the dominant sources of systematic uncertainties in the V+jets measurements is the 
determination of the jet energy, which affects the jet counting. We consider three sources of 
uncertainty related to the jet energy scale: 

• Corrections applied to the measured jet energy to account for the detector response 
and inhomogeneities. The corrections are derived from measured jet px in dijet and 
photon+jet events and are available as a function of rf and px pT|. 

• The dependence of the detector response on the jet flavour. The difference between 
the flavour composition of jets in V+jets events and the flavour composition of the 
jet sample used to extract the corrections is accounted for by an additional 2% un- 
certainty on the jet energy. 

• The average energy removed by the pileup subtraction method includes the activity 
due to the underlying event. We find that the jet energy is systematically decreased 
by 500 MeV in comparison with events without pileup. This amount is included as 
a systematic uncertainty on the jet energy. 

The above uncertainties on the jet energy are added in quadrature and their effect is evaluated 
on the jet multiplicity distribution using simulation. Similar results are found in all the chan- 
nels, for both W and Z events. For W events, the effect of the mismeasured jet energy on Mj is 
evaluated in the fitting procedure. 

Uncertainty on the jet energy resolution also affects the jet multiplicity. The resolution is un- 
derestimated in the simulation by (10 ± 10)% |40|. The effect of this uncertainty on the jet 
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multiplicity has been studied in simulated W + jets events and found to be below 2%. 

The pileup subtraction method was tested by comparing the jet multiplicity in two simulated 
signal samples, one with pileup and one without, with the pileup subtraction applied in the 
former case. The difference, due to residual effects and jets in pileup events, is below 5%. 

All the above uncertainties are added in quadrature and the resultant contribution to the rela- 
tive systematic uncertainty on jet rates is shown in Tables |2] and |3j Also shown in these tables 
are the uncertainties associated with the selection efficiency and the signal extraction proce- 
dure. The largest source of uncertainty related to signal extraction in W events comes from the 
uncertainty of the b-tagging efficiency and mistag rates. This uncertainty is estimated with a 
control sample of tt events decaying into e^f/^ and a pair of b jets. The uncertainty due to the 
selection efficiency in Tables |2] and |3] includes the uncertainty on the lepton efficiency and on 
the selection procedure. As discussed in Section |6| it is evaluated with Z + jets data samples, 
with different strategies for the electron and muon channels. This difference results in a larger 
systematic uncertainty on the muon channel. While the systematic uncertainty on the energy 
scale is correlated among the different jet multiplicities, all other uncertainties are not. The 
relative statistical uncertainty is also shown. 

All statistical and systematic uncertainties are propagated in the unfolding procedure. Uncer- 
tainties due to the unfolding procedure itself are calculated as the differences between the un- 
folded rates using the SVD and the iterative algorithms, and the MadGraph, pythia, and Z2 
or D6T tunes of MadGraph for the unfolding matrix. The resulting uncertainties are shown 
with the final results in the next section. 



10 Results 

Results are given for the leptons within the acceptance defined by the kinematical selection 
cuts in the electron (muon) channels: leading-lepton pj > 20 GeV and \rj\ < 2.5 (2.1); second- 
lepton Pt > 10 GeV and \r]\ < 2.5 (2.4). In addition, electrons in the 1.4442 < < 1.566 
region are excluded, and jets are not counted if AR < 0.3 with respect to an electron from the 
W or the Z. Since the acceptance is different, we do not combine the results for the electron and 
muon channels. 

Two sets of ratios from the unfolded jet multiplicity distributions are calculated The first set 
of ratios is (7(V+ > n jets)/cr(V), where cr(V) is the inclusive cross section, and is presented 
in Tables [i] and [s] and in the upper frames of Figs. [lO- 13 The second set of ratios is cr{Y+ > 



n jets) / (7"(V+ > (n — 1) jets), reported in Tables [6] and |7| and in the lower frames of Figs. 10 - 13 
The contributions of the systematic uncertainties associated with the jet energy scale (including 
pileup effects) and the unfolding are given in the tables and shown as error bands in the figures. 
For n > 2 jets, the pythia pure parton shower simulation fails to describe the data, while 
the MadGraph simulation agrees well with the experimental spectrum. Because of the jet 
threshold Ej > 30 GeV, the sensitivity to the tuning of the underlying event is negligible. The 
expectations from the simulation with the Z2 and the D6T tunes are identical. Thus, in the rest 
of the paper we will only consider the Z2 tune for comparison with the data. 



The statistical uncertainty quoted in the third column of Tables|4]-[7]includes only the statistical 
contribution from the fit results. It is combined (fourth column) with the systematic uncer- 
tainties from the fit and the statistical and systematic uncertainty on the efficiency, which are 
uncorrelated between the samples with different numbers of jets. To estimate the jet energy 
scale uncertainty, jet rates from the fits were scaled higher and lower according to this uncer- 



15 



tainty. Those numbers were then unfolded and the difference in the output from the actual fit 
value is quoted in the fifth column. Finally, the uncertainty due to the unfolding algorithm is 
evaluated as explained above and is shown in the last column of the tables. 



The ratios of W + jets and Z + jets cross sections are shown in Fig. 14 Many important sys- 
tematic uncertainties, such as those on integrated luminosity and jet energy scale, cancel in the 
ratio. The most significant remaining systematic uncertainty is due to the selection efficiency, 
which could be correlated between different bins because of the Mj cut in the W candidate 



Table 2: Relative systematic and statistical uncertainties on the measured jet multiplicity in W 
events, as a function of the jet multiplicity for electron and muon samples. 



Uncertainties on jet rate in W — > ei/ events [%] 


Jet multiplicity 





1 


2 


3 


> 4 


Energy scale and pileup 


±5 


±8 


+11 

-10 


+14 
-12 


+16 
-15 


Selection efficiency 


±0.5 


±0.3 


±1.0 


±1.7 


±4 


Signal extraction 




±0.1 


±0.4 


±3 


±9 


Total systematic uncertainty 


±5 


±8 


+11 

-10 


+14 
-12 


±17 


Statistical uncertainty 


±0.3 


±1.0 


±2.4 


±10 


±28 


Uncertainties on jet rate in W — > /^i/ events [%] 


Jet multiplicity 





1 


2 


3 


> 4 


Energy scale and pileup 


±5 


±8 


+11 
-10 


+14 
-12 


+16 
-15 


Selection efficiency 


±3 


±6 


±4 


±10 


±17 


Signal extraction 




±0.1 


±0.4 


±3 


±9 


Total systematic uncertainty 


±6 


±10 


+13 
-12 


+19 
-17 


±26 


Statistical uncertainty 


±0.3 


±0.9 


±2.4 


±5.6 


±21 



Table 3: Relative systematic and statistical uncertainties on the measured jet multiplicity rates 
in Z events, as a function of the jet multiplicity for electron and muon samples. The signal 
extraction uncertainty is not shown since it is negligible for Z events. 



Uncertainties on jet rate in Z — > e+e events [%] 


Jet multiplicity 





1 


2 


3 


> 4 


Energy scale and pileup 
Selection efficiency 


±5 
±1.1 


±8 
±1.0 


+11 
-10 

±2.2 


+14 
-12 

±2.6 


+16 
-15 

±6 


Total systematic uncertainty 


±5 


±8 


+11 
-10 


+14 
-12 


+17 
-16 


Statistical uncertainty 


±1.1 


±3.1 


±7 


±17 


±43 


Uncertainties on jet rate inZ ^ fi^fi events [%] 


Jet multiplicity 





1 


2 


3 


> 4 


Energy scale and pileup 
Selection efficiency 


±5 
±3 


±8 

+6 

-5 


+11 
-10 
+7 
-6 


+14 
-12 

±10 


+16 
-15 
+24 
-12 


Total systematic uncertainty 


±6 


±10 


+13 
-12 


+18 
-16 


+30 
-21 


Statistical uncertainty 


±0.9 


±2.6 


±5.2 


±18 


±41 
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Figure 10: The ratios cr{W+ > n jets)/(7(W) (top) and cr{W+ > n jets)/cr(W+ > (n - 1) jets) 
(bottom) in the electron channel compared with the expectations from two MadGraph tunes 
and PYTHIA. Points with error bars correspond to the data. The uncertainties due to the energy 
scale and unfolding procedure are shown as yellow and hatched bands, respectively. The error 
bars represent the total uncertainty. 
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Figure 11: The ratio cr(W+ > n jets)/(7(W) (top) and (7(W+ > n jets)/c7-(W+ > (n - 1) jets) 
(bottom) in the muon channel compared with the expectations from two MadGraph tunes 
and PYTHIA. Points with error bars correspond to the data. The uncertainties due to the energy 
scale and unfolding procedure are shown as yellow and hatched bands, respectively. The error 
bars represent the total uncertainty. 
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Figure 12: The ratio a{Z+ > n jets)/c7"(Z) (top) and a{Z+ > n jets)/t7"(Z+ > {n — 1) jets) 
(bottom) in the electron channel compared with the expectations from two MadGraph tunes 
and PYTHIA. Points with error bars correspond to the data. The uncertainties due to the energy 
scale and unfolding procedure are shown as yellow and hatched bands, respectively. The error 
bars represent the total uncertainty. 
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Figure 13: The ratio (t{Z+ > n jets)/(7(Z) (top) and a{Z+ > n jets)/t7"(Z+ > {n — 1) jets) 
(bottom) in the muon channel compared with the expectations from two MadGraph tunes 
and PYTHIA. Points with error bars correspond to the data. The uncertainties due to the energy 
scale and unfolding procedure are shown as yellow and hatched bands, respectively. The error 
bars represent the total uncertainty. 
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selection. The maximal difference observed between the measured and expected values is at 
the level of one standard deviation, neglecting the uncertainties on the theoretical predictions. 
The difference between the expected value of the ratio in the electron and muon channels is 
due to the larger electron acceptance in tj. 



The charge asymmetry, defined as Aw 



t7-(w+)+t7-(w-) ^ measured as a function of the number 
of jets in the events for the muon and electron channels, by fitting separately the events with 
positive and negative lepton charge. The sample size does not allow a measurement in events 
with four jets or more. Table |8] and Fig. 15 show the measured value of the charge asymme- 
try for the electron and muon events. The systematic uncertainties include those from the jet 



Table 4: Results for o'{W+ > n jets)/(7(W) in the electron and muon charmels. A full descrip- 
tion of the uncertainties is given in the text. 



n jets 


a{W+>n jets) 

(7{W) 


Stat. 


Stat. + efficiency 
and signal extraction 


energy 
scale 


unfolding 


electron channel 


> 1 jets 


0.133 


0.002 


0.002 


+0.019 
-0.017 


±0.001 


> 2 jets 


0.026 


0.001 


0.001 


±0.004 


±0.001 


> 3 jets 


0.0032 


0.0004 


0.0004 


+0.0006 
-0.0005 


±0.0001 


> 4 jets 


0.00056 


0.00017 


0.00018 


+0.00012 
-0.00010 


+0.00006 
-0.00001 


muon channel 


> 1 jets 


0.136 


0.002 


0.007 


+0.019 
-0.017 


±0.001 


> 2 jets 


0.026 


0.001 


0.002 


±0.004 


+0.002 
-0.001 


> 3 jets 

> 4 jets 


0.0041 
0.00059 


0.0003 
0.00011 


0.0005 
0.00017 


+0.0008 
-0.0006 
+0.00012 
-0.00010 


+0.0003 
-0.0001 
+0.00001 
-0.00015 



Table 5: Results for (7(Z± > n jets) / c(Z) in the electron and muon channels. A full description 
of the uncertainties is given in the text. 



n jets 


a{Z+>n jets) 
a{Z) 


Stat. 


stat. + efficiency 
and signal extraction 


energy 
scale 


unfolding 


electron channel 


> 1 jets 


0.151 


0.006 


0.006 


+0.021 
-0.019 


±0.001 


> 2 jets 


0.028 


0.003 


0.003 


±0.004 


±0.001 


> 3 jets 

> 4 jets 


0.0039 
0.00070 


0.0009 
0.00036 


0.0009 
0.00036 


+0.0007 
-0.0006 
+0.00014 
-0.00012 


+0.0003 
-0.0001 
+0.00005 
-0.00004 


muon channel 


> 1 jets 


0.149 


0.005 


0.011 


+0.022 
-0.020 


±0.001 


> 2 jets 


0.027 


0.003 


0.004 


±0.004 


±0.001 


> 3 jets 

> 4 jets 


0.0042 
0.00087 


0.0011 
0.00050 


0.0012 
0.00056 


+0.0008 
-0.0006 
+0.00017 
-0.00015 


+0.0001 
-0.0003 
+0.00010 
-0.00001 
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energy scale, the charge misidentification, and the positive versus negative lepton efficiency 
difference. The results are compatible with the inclusive charge asymmetry measurement |28( . 

The charge asymmetry depends on the number of associated jets because the fraction of u (d) 
quarks contributing to the process is different in each case. The measured values are found to be 
in good agreement with the predictions based on MadGraph with the Z2 tune, while pythia 
does not describe well the W charge asymmetry, even for events with a single associated jet. 

Finally, we test for Berends-Giele scaling [13J and measure its parameters. Events are assigned 
to exclusive jet multiplicity bins (inclusive for n > 4), and the corrected yields are fitted with 

Table 6: Results for a{W+ > n jets) /cr(W+ > {n — 1) jets) in the electron and muon channels. 
A full description of the uncertainties is given in the text. 



n jets 


a{W+>n jets) 
a(W+>(n-l) jets) 


Stat. 


stat. + efficiency 
and signal extraction 


energy 
scale 


unfolding 


electron channel 


> 1 / > jets 

> 2 / > 1 jets 

> 3 / > 2 jets 

> 4 / > 3 jets 


0.133 
0.195 
0.125 
0.173 


0.002 
0.007 
0.014 
0.046 


0.002 
0.007 
0.015 
0.049 


+0.019 
-0.017 
+0.002 
-0.001 

±0.004 

+0.003 
-0.004 


±0.001 

+0.012 
-0.001 
+0.002 
-0.004 
+0.017 
-0.003 


muon channel 


> 1 / > jets 

> 2 / > 1 jets 

> 3 / > 2 jets 

> 4 / > 3 jets 


0.136 
0.190 
0.160 
0.144 


0.002 
0.005 
0.011 
0.025 


0.007 
0.014 
0.018 
0.037 


+0.019 
-0.017 
+0.004 
-0.003 
+0.004 
-0.003 
+0.002 
-0.003 


±0.001 

+0.016 
-0.001 
+0.004 
-0.002 
+0.001 
-0.043 



Table 7: Results for a{Z+ > n jets)/(7'(Z± > (n — 1) jets) in the electron and muon channels. 
A full description of the uncertainties is given in the text. 



n jets 


a{Z+>n jets) 
(7(Z+>(f!-l) jets) 


Stat. 


stat. + efficiency 
and signal extraction 


energy 
scale 


unfolding 


electron channel 


> 1 / > jets 

> 2 / > 1 jets 

> 3 / > 2 jets 

> 4 / > 3 jets 


0.151 
0.185 
0.138 
0.181 


0.006 
0.017 
0.030 
0.085 


0.006 
0.017 
0.030 
0.085 


+0.021 
-0.019 
+0.002 
-0.001 

±0.004 

+0.003 
-0.004 


±0.001 

+0.006 
-0.001 
+0.008 
-0.003 
+0.014 
-0.021 


muon channel 


> 1 / > jets 

> 2 / > 1 jets 

> 3 / > 2 jets 

> 4 / > 3 jets 


0.149 
0.180 
0.158 
0.207 


0.005 
0.016 
0.036 
0.104 


0.011 
0.023 
0.043 
0.117 


+0.022 
-0.020 

±0.003 

+0.002 
-0.001 
+0.002 
-0.003 


±0.001 

+0.011 
-0.001 
+0.001 
-0.017 
+0.031 
-0.001 
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Figure 14: Ratio of the W + jets and Z + jets cross sections for the electron channels (left) and 
the muon channels (right) as a function of the jet multiplicity. The ratios are normalized to the 
inclusive W/Z cross section. The expectations from MadGraph and PYTHIA simulations, both 
with the Z2 tune, are shown. These expectations do not differ significantly and they are both 
in agreement with data. The difference between the expected values of the ratio in electron 
and muon channels is due to the larger electron acceptance in t]. Error bars on data are shown 
for the statistical and total uncertainties. The uncertainties due to the energy scale and the 
unfolding procedure are shown as yellow and hatched bands, respectively. 



Table 8: W charge asymmetry Aw in the data and in the MadGraph and pythia (Z2 tune) 
simulations for the electron and muon channel. The uncertainties on the simulation values are 
statistical only. 



n jets 


data 


MadGraph Z2 


PYTHIA Z2 


electron channel 




>o 


0.217 ± 0.004 (stat.) ± 0.006 (syst.) 


0.228 ± 0.001 


0.216 ± 0.003 


> 1 


0.179 ± 0.010 (stat.) ± 0.007 (syst.) 


0.179 ±0.004 


0.267 ± 0.007 


> 2 


0.16 ± 0.03 (stat.) ± 0.01 (syst.) 


0.183 ±0.010 


0.281 ± 0.020 


>3 


0.17 ± 0.08 (stat.) ± 0.06 (syst.) 


0.19 ±0.02 


0.33 ± 0.05 


muon channel 




>o 


0.223 ± 0.003 (stat.) ± 0.010 (syst.) 


0.224 ± 0.001 


0.237 ± 0.003 


> 1 


0.175 ± 0.010 (stat.) ± 0.011 (syst.) 


0.179 ±0.003 


0.222 ± 0.008 


> 2 


0.18 ± 0.03 (stat.) ± 0.02 (syst.) 


0.190 ±0.008 


0.273 ± 0.023 


>3 


0.22 ± 0.07 (stat.) ± 0.02 (syst.) 


0.19 ±0.02 


0.26 ± 0.06 
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Figure 15: W charge asymmetry Aw versus the inclusive jet multiplicity. Data are compared 
with predictions from the MadGraph and PYTHIA simulations. Error bars on the data points 
show the statistical and total uncertainties. Error bars on the simulation points correspond to 
the statistical uncertainty only. Left: electron decay channel, right: muon decay channel. 



the assumption that they conform to a scaling function: 



Cn 



O'n+l 



a, 



(1) 



where (/„ = cr(V+ > n jets), and a is a constant. Previous measurements have shown that the 
Berends-Giele ratio can be approximately constant [7,i9i|41J- Phase-space effects, however, can 
violate this simple proportionality. Therefore we introduce a second parameter, /3, to allow for 
a deviation from a constant scaling law: 



Cm 



(2) 



Because of the different production kinematics of the n = sample, where no reconstructed 
jets recoil against the vector boson, the scaling expressed in Eq. (j2]| is not expected to hold, 
so we do not include the n = sample in the fit. The results of the fit for a. and (5 from the 
W + jets and Z + jets samples are shown in Fig. [l6] The values are given in the (a, j6) plane 
and are compared with expectations from MadGraph with the Z2 tune for the underlying 
event description. Generator results are obtained at the particle jet level for events within 
the lepton and jet acceptance of the selection. The electron and muon expected values differ 
mostly because of the AR > 0.3 requirement between the jets and the leptons, which is applied 
only in the electron channel. The ellipses correspond to 68% confidence level contours using 
the statistical uncertainty only. The arrows show the displacement of the central value when 
varying each indicated parameter by its estimated uncertainty. The fit results are also reported 
in Tables |9] and [lOj The data are in agreement with expectations within one or two standard 
deviations depending on the channel. Furthermore, the Berends-Giele scaling hypothesis i.e., 
the relationship shown in Eq. ([T|, is confirmed to work well up to the production of four jets. 
The /3 parameter lies within one standard deviation from zero for the W + jets case and within 
0.5 standard deviations for the Z + jets case. 



11 Summary 

The rate of jet production in association with a W or Z boson was measured in pp collisions 
at ^/s = 7 TeV. The data were collected with the CMS detector in 2010 and correspond to an 
integrated luminosity of 36 pb^^. The W + jets and Z + jets samples were reconstructed in the 



22 



11 Summary 



CMS 



==^1.5 



T 1 

a(W + > n jets) 



a{W + > (n+1)jets) 



: a + p X n 



uncertainties: 
■ efficiency 
— > jet energy scale 

MG+D6T migration matrix 



0.5 



- 




W ^ ev 

• data 

o MadGraphZ2 'v' 
Stat, only 68% C.L. contours 



36 pb"^ at \[s = 7 TeV 
1 



uncertainties: 
— ► efficiency 
—t> jet energy scale 

MG+D6T migration matrix 




• data 

o MadGraph Z2 
Stat, only 68% C.L. contours 



1.5 



0.5 



■- 




Figure 16: Fit results for the Berends-Giele scaling parameters a. and /3 after pileup subtraction, 
efficiency corrections, and unfolding of detector resolution effects: (top) W + jets, (bottom) 
Z + jets, (left) electrons, (right) muons. The data are compared with the expectations from the 
MadGraph simulation with the Z2 tune. The ellipses correspond to 68% confidence level 
contours considering the statistical uncertainty only, for both data and simulation. The ar- 
rows show the displacement of the central value when varying each indicated parameter by 
its estimated uncertainty. The arrows labelled "MG+D6T migration matrix" correspond to the 
displacement when MadGraph simulation with the D6T tune is used for the unfolding. 
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Table 9: Results for the Berends-Giele parameters in the electron channel compared with ex- 
pectations from the MadGraph Z2 simulation at the particle level. The prediction uncertainty 
is statistical only. 







data 


Stat. 


energy scale 


efficiency 


tune 


theory 


w 
z 


a 

a. 


4.5 
0.6 
4.4 


±0.5 
±0.4 
±0.9 


+0.1 
-0.3 
+0.1 
-0.4 
+0.8 
-0.2 


+0.2 
-0.4 
+0.2 

-0.1 

±0.05 


+0.1 
-0.3 
+0.2 


5.20 ± 0.05 
0.20 ± 0.04 
5.3 ±0.1 






1.1 


±0.8 


+0.6 
-0.4 


±0.3 


+0.1 


0.17 ±0.07 



Table 10: Results for the Berends-Giele parameters in the muon channel compared with expec- 
tations from the MadGraph Z2 simulation at the particle level. The prediction uncertaiaty is 
statistical only. 







data 


Stat. 


energy scale 


efficiency 


tune 


theory 


w 

z 


OL 


4.6 
0.7 
6.4 


±0.3 
±0.3 
±1.2 


+0.1 
-0.3 
+0.3 
-0.2 
+0.1 
-0.3 


+0.2 
-0.1 

±0.3 
±0.1 


-0.5 
+0.3 
-0.3 


5.17 ± 0.09 
0.22 ± 0.07 
4.8 ± 0.1 






-0.5 


±0.9 


+0.1 
-0.2 


±0.2 


+0.1 


0.34 ± 0.09 



electron and muon decay channels. The measurement was performed using particle-flow jets 
with Et > 30 GeV and clustered with the anti-fcx algorithm with a size parameter K = 0.5. 

The Et spectra agree well with the predictions from simulations based on MadGraph inter- 
faced with PYTHIA, and using the Z2 tune for the underlying event description. 

Detector resolution effects were unfolded to extract the exclusive jet multiplicity distributions 
and to measure the ratios of the normalized inclusive rates cr{'V+ > n jets)/c7'(V). The ratio of 
the W + jets to Z + jets cross sections, [(7(W± > n jets)/cr(Z± > n jets)]/ [(7(Z)/(7(W)], and 
the W charge asymmetry. Aw, were also measured as functions of the jet multiplicity. Finally, 
a quantitative test of Berends-Giele scaling, parametrized as a function of two parameters de- 
termined by a fit, was performed. 

All results are in agreement with the predictions of the MadGraph generator, using matrix- 
element calculations for final states with jets, matched with PYTHIA parton shower. In contrast, 
the simulation based on parton showers alone fails to describe the jet rates for more than one 
jet, and the W charge asymmetry, even in the case of only one jet. 
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